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A two-parameter family of quantum spin ladders with local bilinear and bi- 
quadratic interactions is shown to be solvable by a mapping onto fragments of 
integrable spin 1 chains. The phase diagram, consisting of four phases, and the 
ground state properties are discussed. In one novel phase, the ground state is made 
up of plaquette singlets and rung singlets, alternating with a three-rung periodicity. 
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I. INTRODUCTION 

Integrable models have provided extremely valuable insight into the physics of one dimensional 
quantum spin chains. There is a substantial lack of exact results for analogous two dimensional 
systems and the solution of quantum spin ladders might provide a first step in that direction. 
Besides, spin ladders are currently experimentally accessible, interesting in their own and possibly 
related to high Tc superconductivity j^]. 

It has been established that the basic n-leg Heisenberg ladders, with bilinear exchange inter- 
actions along rungs and legs, are gapful (spin- liquid state) for n even and gapless for n odd [|l|. 
These ladders do not seem to be integrable. On the other hand, examples of integrable ladders, 
containing additional biquadratic interactions, have been found and solved by some form of Bethe- 
ansatz (BA) [^,^. Luckily, biquadratic interactions do arise in physically realizable systems, and 
a large class of these generalized, but still 5[/(2)-invariant ladders, have been proven to have a 
matrix-product (MP) ground state Still, the MP approach determines the ground state but, 
with the exception of few lucky cases, the whole set of excitations remains unknown. 

In this paper, the following two-parameter family of 5L/(2)-invariant ladder hamiltonians H = 
Hk=i Hk,k+i will be studied 



*E-mail Giuseppe.Albertini@mi.infn.it 



1 



-2(^(Sfe • Sfe+i)(tfc • tfc+i) + (Sfe • tfe+i)(sA;+i • tk)j + jj.l{Sk ■ tk + Sk+l ■ tfc+i) 

+/i2(Sfc • tfc)(Sfe+i • tfe+i) (1) 

and the eigenvalue problem solved for any /ii, fl2- Here and are spin 1/2 matrices sitting on 
the first and the second leg, respectively. Hamiltonian is actually a special point in a wider 
three-parameter class for which the relative strength of the first two terms is left arbitrarty, (see 
(pi])). All points in this manifold have the property of being reducible, through the introduction 
of the composite rung spin 

Sk = Sk+ tk (2) 

to the general spin 1 chain 

N 

W) = J2 (cos^Sfe • Sk+i+sm9{Sk ■ Sk+i)^) (3) 

fc=i 

In particular (|l|) corresponds to the BA solvable purely biquadratic spin 1 chain, 6 ~ -~n/2. 
Selecting 6 — 0, the Heisenberg spin 1 chain, would remove the second brackets in (|l]) This 
seems physically more appealing but prevents the detailed analysis allowed by BA, since such chain 
is not integrable. Details about definitions and exact diagonalization are the subject of sections II 
and III. 

The system defined by (|l]) has four phases, connected by first order transition lines, namely 
discontinuities in some first order derivative of the ground state energy per site. At least one phase 
is novel because its ground state is a threefold degenerate global singlet where dimerization takes 
place alternatively along legs and rungs with a three-step periodicity. A second phase has a huge 
(^ 3^/2) ground state degeneracy. All this is discussed in section IV. 

In section V, exact elementary excitations are examined in the phase whose ground state is, 
effectively, that of the biquadratic chain, also a global singlet but made up of rung triplets. The 
simplest ones are created by introducing one rung singlet and leaving the remaining iV — 1 triplets 
locked into the ground state of the biquadratic chain with free ends. If iV — 1 is odd, such ground 
state is presumably an SU (2) triplet containing one dynamical kink and one finds a band of 
singlet-triplet excitations, depending on two degrees of freedom, whose energy is degenerate in the 
singlet position. Such excitations had been found in a related ladder but their dependence 
on two degrees of freedom had not been discussed. On the other hand, similar ladders including 
biquadratic interactions are known to have two-parameter singlet-triplet excitations , but their 
nature is different, being a pair of kink-antikink over a dimerized ground state, which is not what 
happens here. In section VI, some general results are discussed for the wider three-parameter 
family of spin ladder hamiltonians. 
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For the sake of completeness, the Appendix contains a review of known results on the biquadratic 
chain. 



II. DEFINITION OF BASIC INVARIANTS 

Each elementary plaqucttc involves four spins, say (si, ti, S2, t2). If only SU{2) invariancc is 
required, the most general hcrmitean plaquettc hamiltonian iJi.2 is a linear combination with real 
coefficients of 14 hermitean invariants. If, beside SU(2) invariance, one requires symmetry under 
the exchange of the two legs and symmetry under the exchange of the two rungs, symmetries 
implemented by operators C and P 

CsiC = ti (i = l,2) C2 = l 

PsiP = S2 PtiP = t2 P^ = 1 

then Hi^2 is a linear combination with real coefficients of 8 hermitean invariants. 

The proof goes as follows. The local plaquette Hilbert space is 
C^^. It breaks into the orthogonal sum of SU{2) multiplets: one quintuplet (dim=5), three triplets 
(dim=9) and two singlets (dim=2). Within each multiplet, all states are obtained by application 
of the lowering operator the highest weight vector (h.w.v.) v^^\ where 

Si,2=Si+S2+ti+t2 S±2 = ST,2 ± *Sf ^2 

S+2t;(^) = Sl^v'-"^ = s{s + v'-"') e Hi,2 

ill, 2 is 5t/'(2)-invariant, i.e. [iii,2, 81,2] = 0, therefore it maps h.w.v. into h.w.v. of the same spin. 
So, spin-2 h.w.v. must be mapped into itself, each spin-1 h.w.v. will in general be mapped into 
a linear combination of all three spin-1 h.w.v. and the same happens for the two spin-0 h.w.v.. 
Altogether, imposing hermiticity, 6 real parameters plus 4 complex (hence 8 real) ones. The total 
spin Si, 2 commutes with C and P, so h.w.v. can be labelled by C and P eigenvalues. Define, with 
self-explanatory notation 

eo = |TT>« /o = |TT>t eo = Ui>. 7o-Ui>t 
ei = UT>« /i = UT>t e2 = |n>. /2 = |n>t 

The 6 h.w.v. v^^^ will be chosen to be 

v^^^ = eofo (2,1,1) 

^^1^^ = eo/i - eo/2 + ei/o - 62/0 (1, 1, -1) 

4^' = eo/i + eo/2 - ei/o - 62/0 (1, -1, 1) 

= eo/i - eo/2 - ei/o + 62/0 (1,-1,-1) 
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4°^ = ei/i + 62/2 - ei/2 - e2/i (0,1,1) 



.,(0) - 



ei/i + ei/2 + 62/1 + 62/2 - 2eo/o " 2eo/o (0, 1, 1) 



The three numbers on the right column are eigenvalues of S'1,2, C , P. If C and P symmetries are 
imposed. Hi 2 maps v^'^^ into itself and each wj^-* {i — 1,2,3) into itself, while its action on the 
span of {vf''' ] i = 1, 2} is determined by 4 real numbers 

iJi^2i^(2) = to(2)^;(2) ^ mp)i;f ) z = 1, 2, 3 
2 

This proves that the general i?i,2 is a linear combination with real coefRcients of 8 linearly 

independent invariants (notice incidentally that h.w.v. v^^^ are orthogonal but unnormalized, 
Ij-yj*^'*!!^ — so rri{2 — 3"i2i''*). One can choose 7 of them to be three bilinear (Heisenberg) 
terms 

I^ll = Si • S2 + ti • t2 (4) 

I'^ll = Si • ti + S2 • t2 (5) 

if^ = Si • t2 + S2 • ti (6) 

three biquadratic (plaquette) terms 

/j^2^=(si.S2)(ti.t2) (7) 

/j;^2^ =(si.ti)(s2-t2) (8) 

=(si.t2)(s2-ti) (9) 

and the identity I. An eighth one is necessary to have a complete set, but it involves more 
complicated combinations of the basic spins and it will not be needed in the following. To show 
that the six in (|4|)-(^ plus the identity are indeed linearly independent, write 

6 

,2 



i7i,2 = col + ^ Ckl^^^ 



k=l 



It is a matter of easy algebra to find the action of the basic invariants (0)-(^) on the h.w.v., resulting 



I2-] . Ci+ C2+C3 , C4 + C5 + Cg 

' = co-\ — 



m\' = Co 



2 16 

Cl + C2 — C3 — 3C4 + C5 — 3C6 



2 16 



(1) , Cl - C2 - C3 C4 - 3C5 - 3C6 
-2 =^^0 + ^ + ^ 

= + -C1-C2 + C, ^ -3c4-3c, + ca ^^^^ 
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(0) — 3ci 9c4 + 3c5 + 3cq 
mn =^'0 + ^ + 

(0) , Ci - 2C2 - 2C3 C4 + 7C5 + 7C6 

m^2 = Co + + 
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(0) _ C2 - C3 -C5 + C6 
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and = Stoj'^''. The vanishing of aU {m'f^} implies the vanishing of aU {q}, proving hnear 
independence of the invariants chosen. 

Scalar products of the rung spin will now be expressed in terms of these invariants. Clearly 

<S q - 7-(l) , r(3) 
Si • &2 — -'1,2 + ^1,2 

while the biquadratic term (Si • 82)^ seems to involve more complicated invariants. But since it 
is invariant under C and P, it must be a linear combination of the 8 basic ones. Actually, the six 
(@)"@ the identity are sufficient: from ( [l0| ) 

(Si . 82)^ = f - ^ + 1^2^ - ^ + 21^ + 2li$ 

Consider now the one-parameter spin ladder Hamiltonian (|^) where each spin S is the composite 
object defined in (^). Two parameters can be added after noticing that the charges 

N N 

Qi = ^^^k Q2 = SfcSfc^i 

k=l k=l 

commute with Hq{9) for any 9. This can easily be checked by a direct calculation. Finally, the 
Hamiltonian to be studied is 

N 



E 

k=l 



sm 



(cose- + 4'^+i) + 2sin^^(45+l + + {sme + Ml + 3^2)45+1 



+4M2/g+i 



+ A^(^ + ^ + ^) (11) 
4 2 4 



Eq. (|l|) is dill), up to a constant shift, if one takes = — 7r/2 and 

Ml = -1 + Ml + 3^2 A2=4m2 (12) 

In the following, parameters (mi,M2) will be adopted to describe the model. It is always possible, 
through (p^, to revert to the original ones (mi,M2)- 

One further remark. Since each rung space V^^"' V^*^ carries a spin 1 ® representation, not 
just spin 1, (Sfc • Sfc+i)'^ is hnearly independent from Sk ■ Sk+i and (Sfc • Sfc+i)^. But, from (p^, it 
can be seen that (Sk ■ Sfc+i)'^ is a linear combination of the first two powers, X, and Ik^l^i-, 

that is why adding a cubic term to Ho{9) does not generalize ([ll|). 
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The composite rung spin has, of course, been introduced in several previous works. It lies behind 
the physical idea that even-legged and odd-legged Heisenberg ladders should be in different phases 
101 . More closely to this work, it has been used Q to map two-leg ladders into the AKLT chain |l3) 
((||) at tan 6' = 1/3). Earlier, Xian and then Kitakani and Oguchi ^] had studied ( pi] ) at 6* = 0, 
but without the rung-rung interaction ij^^l^i which is actually responsible for the rise of two new 
phases. 

III. FRAGMENTATION AND BETHE-ANSATZ 

Charges Qi and Q2 have a simple physical meaning. In each four dimensional rung space 
Vk ~ V^'"'' ® ^fe*' — OJ^e introduces the singlet-triplet basis {\s >k] \t >k,t = —1,0, 1} consisting 
of the singlet and the triplet of the rung spin (||), i.e. 

Sl\s>k=0 Sl\t>k=2\t>k S^\t>k^t\t>k 

Then ^ counts the number of triplets and ^ counts the number of pairs of neighboring triplets. 

Now consider the 3^ dimensional subspace (of the total 4^- dimensional Hilbert space of the 
ladder) spanned by |ti,t2, ....,ijv >i {tk = ^1,0, 1). In this sector Qi and Q2 are constant, at 2N 
and AN respectively, and Hq acts effectively as a spin-1 chain with periodic boundary conditions 
(p.b.c). Such chain is Bethe-ansatz solvable in three cases: (a) 9 — j, (b) 6 — ~j and (c) 
6* = — -I In case (a) it is the S'C/(3)-invariant Sutherland-Uimin chain |jri|, in case (6) the Babujian- 
Takthajian chain jl^ and in case (c) the purely biquadratic chain [ p^ . 

Next, consider states containing one singlet and — 1 triplets. Eigenvalues of the two conserved 
charges are fixed at 

Qi=2(iV-l) g2=4(A^-2) 

regardless of triplet's position. If the singlet is, say, on the N*^ rung (due to p.b.c. the singlet's 
position is actually immaterial), Hq acts on these vectors like 

N-2 

Ho^Y. (cos^Sfc • Sk+i + sme{Sk ■ Sfc+i)^) (13) 
fc=i 

that is exactly like a spin-1 chain of length — 1 and free boundary conditions (f.b.c). The 
singlet can be positioned anywhere, it just opens a fracture in a ring of spins, consequently each 
eigenvalue of (|l^) appears with an A^-fold degeneracy in the ladder spectrum. More importantly, 
the spectrum of (|l^) can be found exactly by Bethe-ansatz, at least in cases (6) and (c) which 
have been more or less extensively studied for f.b.c. [l^Jl^ . The opposite situation arises when all 
rungs are singlets. There is only one such state and its eigenvalue is trivially zero. 
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Between these extreme cases, and orthogonal to them, there is all possible intermediate situ- 
ations, characterized by alternating fragments of triplets and singlets. Sectors are labeled by a 
sequence of positive integers, the lengths of fragments, {n[^\ n[^\ N2^\ N2^\ Nn \ Nn^}, with 
j^n^^ j^{0) ^ ^(0)^ ^j^^ ^ ^(1) ^j^j ^(0) _^ ^(1) ^ ^_ ^g^^i^ ggj^^Qj. jg spanned by the 3^'" 

vectors 



^ — \ ^: ■■■7 I'l : ''2 ' • •• ' )v(l) ; '^1 ■ 1 '^^ J ''1 I ''2 ''■'' ^(l) 



> 



On all these basis vectors 

n 



)|s,S, ...,S > (0) 



where i/p"''' acts on the string of triplets like a spin-1 chain of length Nj^'^ and f.b.c. Call lip > 
any of its eigenvectors and Ej{Nj^^) the relevant eigenvalue. Then 



\lp >= \S, S, S >j^(0) (^l-ipl >^(l) (8)|S, S, S >^_(0) (g)... (g) IV'n >j^(i) (14) 

is an eigenvector of Hq with eigenvalue 

n 

E = Y,E,iN^}^) (15) 
i=i 

Vectors in ([T^ ) provide a complete set for the whole ladder to the extent the vectors \ipj >Nj 
provide a complete set of eigenvectors for the spin-1 chain of length Nj, so diagonalization of the 
spin 1 chain with p.b.c. and f.b.c. provides a complete solution to the diagonalization problem of 
the spin ladder (pi]). In general this is only a partial simplification. Instead, in cases (a), (6) and 
(c) eigenvalues can in principle be found, for all fragments, by the suitable Bethe-ansatz. What 
happens here is very similar, actually almost identical, to what happens in the mixed Heisenberg 
chains studied by Niggemann et al. Q. There, composite rung spins alternate with single spins 
and, very much like it happens here, fragmentation takes place when one rung spin is in a singlet 
state. In the following I will partly adopt notation and methods introduced in their work. 

IV. THE GROUND STATE PROBLEM 



To identify the ladder ground state one must minimize (^_5|). An easy first step is to choose, 
in (p^), the lowest eigenvalue for each fragment. To fix the notation, denote with E^\N] 1) the 
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lowest eigenvalue of the iV-site spin 1 chain with p.b.c. and with Eq'^\n'; 1) the lowest one for the 
iV'-site chain with f.b.c. Then, in a sector {N^^\ Nj'^'' }^^^ the lowest eigenvalue is 

f E;=i(i?^^''(<^l) + (2Mi+4M2)ivf^-4/.2) 7V«^0,iV 
^ = \ S*^' + (2^*1 + i^l2)N 7V(i) = N (16) 

[ n'-^^ ^ 

Following Q, a spin rung is joined to, say, the right of each triplets fragment, bringing its length 
to J + 1. Set Hj to be the number of triplet fragments of length j + I (i.e. j triplets plus one 
singlet at the edge). Consider, for the time being, only fragmented configurations, i. e. iV^°) > 0, 
specified by {uq, ni, ^2, . . .} which must satisfy 

Af-l 

^ (j + l)n, - N 

Their energy is for rig = N, otherwise 

E{no,ni,...,nN) = 1) + J(2mi + 4/i2) - 4^2)%" 

( f) 

where Eq (1; 1) = by definition. In the limit — > 00 the densities wj = nj/N must be chosen 
to minimize 

+00 



lim — 

N-*oo N 



eiwo,wi,W2, ■■■)=J2 {^o\j-^ 1) + i(2/^i + 4Ai2) - 4^2) (17) 



+ 00 



E(J + 1H^1 (18) 
j=o 

As observed in ||^, the problem (|l7|), (^8|) has arisen in classical statistical mechanics of systems 
with competing interactions ||20[] . Owing to the linearity of ( |l7|) in the parameters {wj}, the 
extremum occurs when one Wj 7^ and all others are zero. Then 

Wj = — ^ Wk = (fc 7^ j) 

and the ground state energy per site is 

^ ^ ^(^.^ ^ i^^^^(j;l)+j(2/ii+4M2)-4A^2 ^^g^ 

The relevant phase is denoted < j >. One is left with the task of determining the minimum in 
the sequence of real quantities (p^), j > where it is understood that e(0) = 0. The minimum 
might be reached at j 00, or e(oo) = eo + 2/_ti + 4/i2 where cq is the ground state energy per site 
of the spin chain, independent from boundary conditions. The thermodynamical analysis outlined 
above does not allow to distinguish between the periodically closed sector N^^'' = N and the open 
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sector where, for instance, one singlet is introduced. In appendix B it is shown in detail that, when 
infj>o e(j) = hmj^oo then the ground state indeed belongs to the sector N'-^'> — N. 

All said so far is true for any value of 9 in (pi]). Yet, a comparison of e{j) requires the knowledge 
of E'^^\j; 1), in principle for any j. The BA solvable cases (a), (&), (c) are perhaps not the most 
physically relevant, but they allow an efhcient computation of E^^\j\ 1). In the following I will 
mostly concentrate on case (c), the purely biquadratic spin 1 chain. The information available for 
this chain is summarized in Appendix A. The relevant results on finite size ground state energies 
are gathered in Table |. Some considerations on models with arbitrary 9 are postponed to section 
VI. 

It can be shown that the (/ii , ^2) plane is divided into four regions, corresponding to four different 
ground state energies per site, hence four different phases. Since and ^2 often appear in the 
combination 2/ii + 4/i2 define 

/Z3 = 2/^1 + 4/X2 

First, seek the (/i2,/i3) values for which e(0) is lowest, that is < e(j), j > 1. It is convenient 
to rewrite this condition adding and subtracting a term containing cq (whose numerical value is 



reported in (Al)) 



4Ai2 < so(j) + j(eo + ^l3) J > 1 so(j) ^^o^^ (j; 1) - jeo 

The sequence so(j) is bounded, so a first necessary condition is cq + /ia > 0. Furthermore, from 
Table 2, so{j) is increasing in the range 2 < j < 51, with a mininmm 

so(2) = E^g^\2; 1) - 2eo = -4 - 2e„ ~ 1.5937 (20) 

What can be said for j > 51? From Appendix A, so(j) has different limits at j — )■ 00 for j even 
(£(5,+)^ and j odd (e'*'^^) but they are both larger than —4 — 2eo. It seems, from Table ||, that 
so(j) at j ~ 51 is already rather close to the asymptotic regime, so one can safely, although not 
completely rigorously, conclude that sq^j), j > 2, is minimal at j ^ 2 and that, since eo + /is > 



inf 

J>2 



(soij) + j(eo + Ma)) = so{2) + 2(eo + Ma) = -4 + 2^3 



Hence, the three necessary and sufficient conditions for being in phase < > are 

eo + /i3 > 

2/i3 - 4^2 > 4 (phase < >) (21) 

^J■3 - 4^2 > 

The third condition comes from e(0) < e(l), considered separately. 

Phase < 1 >. The conditions for its existence, e(l) < e(0) and e(l) < e{j), j > 2 are rewritten, 
from (lit 
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s^U)'^^ -2E(f\j;l)/{j-l) (j>2) 

where si{j) is listed in Table up to j = 51. It is decreasing in this range and it certainly is at 
large j. In fact, from the known numerical values (|Al|), ( A_8) and ( A9), 2eo + 2e^^'^^ < 0; moreover 

SlU) = f + oil/j) = -2eo - ^^^^ + oil/j) ij oo) 

Hence one can conclude that supj>2 si{j) = si(2) = 8 and phase < 1 > appears when 

Ma - 4^2 < ^3 + 4^2 > 8 (phase < 1 >) (22) 

Phase < 2 >. The conditions for its existence are e(2) < e(0), e(2) < e(l) and e(2) < e(j), j > 3, 
or, remembering that Eq'^\2; 1) = —4 

Ai3 - 2^2 < 2 + 4/i2 < 8 Ai3 + 4Ai2 > S2(j) (j > 3) 

S2(j) =^ (34^'(j; 1) + 4(j + l))/(-j + 2) (j > 3) 

The sequence S2(j) is listed in Table || and it is increasing up to j = 51. It certainly is for large j 
where 

3eoj + 3e(--±) + 4{j + 1) , _6eo - 3e(-'±) - 12 , 
S2U = -— t; + o 1/j = -3eo - 4 + ■ + o 1/j) 

because — 6eo — 36^*'^-' — 12 < 0. This time, it is not so straightforward to guess the upper limit 
of 52 (j)- If S2(j) keeps increasing for j > 51 up to the asymptotic regime where one can rely on 
the previous equation, then sup^^g S2{j) — limj^+oo S2(j) — — 3eo — 4 ~ 4.3906 and phase < 2 > 
appears for 

^3-2^2<2 - 3eo - 4 < /i3 + 4/^2 < 8 (phase < 2 >) (23) 

On the other hand, one cannot rule out the possibility that supj>3S2(j) lie slightly above and 
this would leave a tiny window for an additional phase. However, the existence of phase < 2 > is 
unquestionable . 

Finally, phase < oo >, whose ground state is that of the periodic spin 1 chain, shows up when, 

def 

for each j, e{j) < e(oo) — limj^oo e(j). Since e(oo) — eo + fj.3, this means 

eo + ^J-3<0 eo + At3 < (M3 - 4a*2)/2 

eo + Ai3 < so(j) - 4^2 (j > 2) 

We already know that so(j) is minimal at j = 2, so the third condition amounts to eo + fJ-3 + 4/i2 < 
So (2) = —4 — 2eo which, by itself, implies the second inequality. Hence 
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eo + < ^3 + 4^2 < — 4 — 3eo (phase < oo >) 



(24) 



are the necessary and sufficient conditions that guarantee the appearance of phase < oo >. Ap- 
pendix B completes the proof: the ground state is that of the periodic spin 1 chain because, under 



conditions (24), local variations around the sector A^*^^' ^ N only lead to an increase in energy. 

Under the foregoing assumption on the upper limit of sequence S2(j) there is no room for other 
phases because inequalities (^ll)-(^4|), with their separating lines, fill the whole (/^2,M3) plane. The 
ground state energy per site is e(0) — (phase < >), e(l) = {^13 — 4/X2)/2 (phase < 1 >), 
e(2) ~ (—4 + 2/.t3 — 4/i2)/3 (phase < 2 >) and e(oo) — eq + fj,^ (phase < 00 >). It shows first 
derivative discontinuities at the boundaries. 

Here are some general features of the four phases. Define the total spin 

N 

Stot = Sfc Sjoj — 5*10! ("Stot + 1) 

k=l 

In phase < 00 > the ground state is simply that of the periodic biquadratic spin 1 chain. Supposing 
N even, it is a global singlet and, it was conjectured in ||l5|, its degeneracy should be 2. Phase < > 
has certainly a unique ground state and of course S'tot = 0. As to phase < 2 >, the ground state is, 
again, a tensor product of local singlets, hence a global singlet, 5tot — 0. In fact, it is represented 
by the sequence | . . . , s, t, i, s, i, s,t,t,s, . . . > where neighboring triplets \tk,tk+i > are locked 
into the singlet ground state of — (S^ • S^+i)^. Supposing A^ = (mod 3) to avoid problems 
with AFM seams, the ground state is threefold degenerate: rung singlets sit on rungs labelled 
k (mod 3), with the freedom of choice k = 0,1,2. It clearly spontaneously breaks translational 
invariance. Finally, phase < 1 >, described by the sequence | . . . , i, s, t, s, i, s, . . . > is rather 
different because nothing fixes the state of each isolated triplet resulting in a huge degeneracy, 
2 • 3^/^. Not only is translational invariance spontaneously broken, but one is dealing with A^/2 
effectively noninteracting triplets that can combine in a "ferrimagnetic" state S'tot — N/2, or 
combine in local pairs yielding S'tot = (these are only two of the many possibilities). In other 
words, the ground state is not necessarily an eigenstate of Stot- 



V. EXCITATIONS 

The four phases display a rich variety of excitations. Most of them, especially in the fragmented 
phases < >, < 1 > and < 2 > can easily be determined resorting to (|l^), therefore only phase 
< 00 > will be extensively treated in the following. Appendix B gives a complete discussion of 
how, in phase < 00 >, "local" perturbations of the biquadratic ground state can only lead to an 
increase of energy. 

A first kind of excitations, within the sector A^^^^ = A^, are those of the biquadratic chain 
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itself. Supposing N even, they exsist in even number and have a gap (see Appendix A for a more 
exhaustive discussion) 

A£;gap ~ 0.173179 

A second kind of excitation, called TSST (triplet-to-singlet spin flip) in Q, is obtained by 
introducing one singlet, sector N'^^^ — N — 1. There is N ways to do so and the eigenvalue is 
A^-fold degenerate. This implies that linear superpositions of such states produce eigenstates of 
the shift operator, all with the same energy. As a matter of fact, dispersionless excitations have 
already been found in akin ladders j^. The energy difference is, see ( |l6|) 

AE{N) = E^/\n - 1; 1) - E^''\N; 1) - - 4.^,2 
AE = lim AE{N) = eS"-^ - cq - ^3 - 4/i2 

A^— >oo 

From conditions ( p^ it follows that AE > 2eo + 4 + e'^'''^^ ~ 0.1479 confirming that the energy can 
only be increased by fragmentation. At the same time, that numerical value shows that, at least 
near the boundary lines, the energy scale is comparable, actually slightly smaller, than ASgap. 
Therefore fragmentation may be relevant to the low energy physics of the ladder. Despite what 
might seem, such excitations depend on two degrees of freedom. In fact, as discussed in Appendix 
A, the ground state of the fragment of triplets, i.e. the ground state of the open biquadratic chain 
with — 1 (odd) sites, is only the lowest in a continuous band. This adds a second degree of 
freedom to the singlet's position. Therefore, if n is the singlet position 

AE{a, n) = e{a) - e(7r) + e^'"'^^ - eo - /is - 4/i2 a € (0, tt) 

Here e(a) — e(7r), always positive, measures the excitation energy of the single kink in the open 
spin chain and vanishes at a — > tt, see ( |A6| ). The n-dependence is, of course, trivial. Nonetheless, 
it is at least plausible that under a small variation of coupling constants in (|l|), for instance a 
different relative weight between /^^^ and /^^•', or 7*^'*^ and l'-^-' interactions, AE might acquire a 
genuine two-parameter form. Since the triplets fragment is presumably locked into a spin 1 state 
(see Appendix A), Stot = 1 for such states. Triplet-singlet excitations have been found in similar 
ladders, but they arise through a different mechanism (t]]^. 

When two rung singlets are introduced (A^^^^ = A^ — 2), their energy becomes position dependent. 
If they sit on rungs ni, rt2 and d — n2 ~ ni > 1, then from jl^ ) 

AEiN; ni,n2) = E'^/\d - 1; 1) + E^^\n - d - 1; 1) - eI^\n- 1) - 2^i3 - 8fi2 + 4^2(5^4 

If A^ — > 00, keeping ni, 112 fixed {p{d) is the parity of d) 

1) - 2eo -I- e^^'f (''-I)) - 2/i3 - 8/12 {d>2) 
AE{d) = { -3eo + e(^'-) - 2^3 - 8^2 {d - 2) (25) 

= _ 2fl3 - Afl2 (d = 1) 
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Notice that all these excitation energies are exact. It is easy to see that AE{d) > always. 
It represents a sort of "static interaction energy" between singlets, whose distance dependence is 
governed by so{d). The open chain surface energy is different for even and odd lengths, so so{d) has 
a uniform d-dependence only if the parity of d is fixed. When d is even, AE{d) uniformly decreases 
(see Table |^ from a maximum AE{2) to a minimum lim^^oo AE{d) — — 2eo + 2e('*'^) — 2/i3 — 8/i2. 
When d is odd, the behavior is certainly uniform for d > 3, because AE{d) increases from AE{3) 
to a maximum limd^oo AE{d) = 2e^'*'+^ — 2eo — 2/i3 — 8/12- Still, AE{1) is not necessarily lower 
or higher than AE(3), since 

AE{3) - AE{1) = -4 - 2eo - 4^2 



whose sign is left undetermined from conditions (24). The short range behavior can be either 
repulsive or attractive. 

A third possible kind of excitations, named n-TSSF in Q], is produced by the insertion of a block 
of m singlets, represented by \ . . . ,t,t,t, s, s, s, s, i, t,t,t, . . . >. Now, from (p^, supposing N even 
as always 

AE{N; to) = E^^^ (TV - m; 1) - TO/xg - 4/X2 - E^^^^^ {N; 1) 

AE{m) ^ lim AE{N; to) = -(m - l)(eo + Ma) + ef'^P^"" - eo - /^s - 4^2 

which is positive under conditions (p4|). At eo+/i3 = these perturbations determine the instability 
versus phase < >. 



VI. OTHER VALUES OF 6 



It is easy to see that phases < >, < 1 > and < 00 > will be present no matter what value 9 
takes in (pT|). In this general situation, ( p^ will formally be the same but now withi;^-^^(A^j^^;l;e') 
and Elf\N; l;0), the energies of the spin 1 chain (^) with arbitrary 9. The proof of section 4 goes 
through unchanged. One has to find the minimum of 

The conditions for phase < > amount to 

4^2 < soU) + j{eoi9) + A^s) (j > 1) so(j) E^/\j; 1; 9) ~ eo{9)j (26) 

Notice that so(j) is bounded. Set 

so(ji) = inf so(j) 
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Of course ji might be +00. In this case, since Muij^oc so{j) = e^'^'^\ the two possible surface 
energies, one has so(ji) = min{e^'''+\ e^^^"^}. At any rate, so(ji) is finite. A necessary condition 
required by (^6|) is 60(6*) + /is > 0, but since one does not know the exact value ji one can only 
give sufficient conditions for the existence of < >, namely 

eo{0) + fi3>0 (phase < >) 

4/i2 < inf sa{j) + inf i{en{0) + n^) = so(ji) + eo(0) + ^3 

which can always be simultaneously satisfied by suitable values of (/i2,/i3)- 

Likewise, for phase < 00 >, the condition e{j) > eo{d)+fi3, (j > 0), translates into eo(6')+/i3 < 
and eo(0) + /is + 4/t2 < so(j): (j > l)i which can both be fulfilled taking 

eo{&) + H3 <0 eo{9) + A<3 + V2 < so(ji) (phase < 00 >) 

These conditions are necessary and sufficient, and are compatible. 
Finally, the inequalities for < 1 > are rewitten 

^3_4,.2<o ^, + ^^,>. ^^o\rA;0) ^ej ^.^2 

J - 1 

Sequence si(j) is also bounded. Set si{j2) = supj>2 The ensuing necessary and sufficient 

conditions are 

fj.3 -4^2 <0 M3 + 4/i2 > si(i2) (phase < 1 >) 

On the contrary, it is in general impossible to draw conclusions on the existence of other phases. 
For example, phase < 2 > would be present if it were possible to fulfill simultaneously 

4/i2-2M3>4^^(2;l;0) /i3 + 4/i2 < -24^^(2; 1; 0) M3 + 4m2>S2(j) (27) 

S2U) ((j + 1)4^' (2; 1;^) - 34^' (j; l;0))/(j - 2) (j > 3) 

Notice that S2{j) is also bounded. If we set 52(^3) ~ supj>3S2(i), the second and the third 
inequality in ( p7| ) are compatible if and only if 

S2{]2) < -2Eif\2-l-e) (28) 

At least for 6* = 0, 69(6*) has been determined to a great accuracy and one can attempt to estimate 
whether ( p8| ) holds. It is known that 

eo{e = 0) ~ -1.401484 

On the other hand, not many data seem to have been published for Eq {j; 1; 0) even at small j. 
Diagonalization of (^, with free boundaries, can easily be carried out numerically up to 7 sites (as 
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a test, I diagonalized the periodic chain as well and compared the results with those published in 
p3[). It is sufficient to examine the sector S'f^j ~ because, owing to 5f7(2)-invariance, all possible 
eigenvalues are bound to appear there. Now, Eq^\2] 1;0) — —2, while the sequence {s2(j)}J=3 is 
found to have the approximate values {1, 1.9606, 1.8302, 2.0277, 1.9807}. Furthermore 

lim S2{j) = -2 - 3eo ~ 2.20446 

It is of course impossible to draw a rigorous conclusion from these results, but it looks extremely 
likely that inequality (pSl) is true and therefore phase < 2 > is present also for 6* = 0. 



VII. DISCUSSION 

It is useful to compare the ladders considered here with those studied in previous papers. To 
shorten notation, I will denote by |0 ><x> the ground state of phase < x >, or one of them if it 
is degenerate. In the work of Niggeman et al. similar methods were employed but for different 
systems, namely the mixed Heisenberg chains. Perhaps the most closely related ladder is that of 
Xian [Q, which is (|Tl|) for = and /12 ~ 0. Perturbations of this ladder where considered in 
Xian ffiids only two phases, < > and < cx3 > in the present notation, but, since 6 = 0, 
his ground state |0 ><oo> is that of the spin 1 Heisenberg chain. The absence of phases < 1 > 
and < 2 > can be understood when one notices that even for the case considered in this paper, 
9 = — 7r/2, phases < 1 > and < 2 > would be missing without the rung-rung interaction J^^"* 



can be seen by setting /i2 = in (22) and (|23|). 

The ground state of a wide class of 5t/(2)-invariant ladders has the matrix-product (MP) form 
1^ . Such states can be translationally invariant or dimerized 

N N/2 
IV'o"''V) > = Tr Jl |?/'o'^™^(wi>'"2) Tr J]^ 52n-l(ui)52n(u2) 

n— 1 n— 1 

u|s >„ +|0 >„ -\/2|l >„ 
V2|-l>„ ■u|s>„-|0>. 



9n{u) = 



For models with translationally invariant MP ground states, u = and u = 00 are the only 
points which have a significant overlapping with the ladders presented here. At u = the MP 
state is made up of rung triplets and is, effectively, the ground state of the AKLT spin 1 chain. 
That is what one gets for ( pi] ) by setting tan6' = 1/3 and keeping /ii, /i2 small enough. Now, 
|0 >oo is also made up of triplets, but it is that of the biquadratic spin 1 chain, hence a different 
kind. Furthermore, |0 ><o>, |0 ><i> and |0 ><2> all contain rung singlets and are certainly not 
|^(inv)^Q-j ^ Conversely, after suitable normalization, the limit |i/'q"^^(oo) > is made up entirely of 
rung singlets. There is only one such state so it must coincide with |0 ><o>- Of course excitations 
may differ according to the detailed form of the hamiltonian. At any rate, |0 ><i> and |0 ><2> 
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break translational invariance so they can never be of the form \^pQ^'{u) > for any u. For the 
same reason, |0 ><2>, which is invariant under a three-rung translation, can never be of the form 
\i^a^^'^\ui,U2) > either, since the latter is invariant under a two-rung translation. But, within 
the large lowest energy eigenspace of phase < 1 >, at least one vector is in MP form. To show 
it, notice that at finite, non-vanishing values of ui, M2, |^o"^™^(wi, U2) > is a linear combination 
necessarily containig kets with all triplets and the ket with all singlets. The only way to reproduce 
|0 ><i> is to set ui = and to take M2 — s- 00 after proper normalization of g2n{u2). The outcome 
is a linear combination of kets which have singlets on even rungs and triplets on odd rungs, exactly 
like |0 ><!>. It should be observed though the the MP approach only in few cases allow to find 
the whole spectrum, as it is instead possible from BA for hamiltonian (pj). 
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APPENDIX A: THE BIQUADRATIC SPIN 1 CHAIN 

The purely biquadratic spin 1 chain is integrable for periodic and free [p^ , p^ boundary 
conditions. The definition of the relevant elliptic parameters differs in the two cases, as treated in 
the literature. For the sake of readability, a Landen transformation has been applied to the 
results of [p^, in order to make notations compatible. 

Define an elliptic modulus from the relation 

K'{k) 1 3-x/5 

where K{k) and K'{k) arc the elliptic integrals of first and second kind. The ground state energy, 
which is the same in the two cases, is then 

fl + 4y , ) - -2.796863... (Al) 

n—l 



Co 



For the periodic case, energy and momentum of the excitations are given, in the limit N cx), by 

2iy 2v 



i=l 



~e(p) = —K{k)^/k'^ + k'^ sinp k'^ + k^ ^ 1 (A2) 

TT 

There is a gap in the spectrum 

-t-00 ^ 

AEg^p = 2^(0) = 2^K{k)k' = V5 [| ( j-prr)^ - 0-173178... (A3) 
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The free boundary case has been solved by noticing that the spectrum is the same, up to degen- 
eracies, as that of the XX Z spin 1/2 chain with a boundary field |jl^ 

1 ^"^ sinh 
Hxxz = - 2 II (f^fe^fc+i + ^fe^fe+i - cosh7CT|a|+i) + ^^(aj ~ a^) (A4) 
fe=i 

Ex,ci^Exxz-\{N -I) q = e-' 

The mapping, valid for any N, works because both hamiltonians can be written as sums over 



generators of the same Temperley-Lieb algebra ||1^. Consequently, ground state energies of (A4) 
are for any N identical, up to a shift, to ground state energies of the for biquadratic chain with 
free boundaries. In turn, (A4) has been solved by coordinate Bethe-ansatz [Q, so eigenvalues are 
found from solutions of a set of coupled trascendental equations 



1 1 " / 



Q{a; x) —i In 



sinh(a; + f 



N 



sinh(a; — ^ 

E=-{N-1) cosh 7 - 2 sinh 7^ e' (a,; J) 



= 2 arctan(tan — coth x) (A5) 



The function Q{a;x) is defined to be continuous for real a. The {Ij} in ( [A5| ) are positive integers 
0,|8|. The ground state belongs to the sector n = N/2 {N even) or n = (iV - l)/2 {N odd). 
In both cases, in (|A5|), Ij — j, j — 1,2 . . .n. Data in Table || have been obtained by mmicrically 
solving (A5) with this choice of integers. 

Excitations are, in general, gapful (no dispersion relation can be written here because linear 
momentum is not conserved). 

AE = Y €{a^) e{a) = 2 sinh7:^dn(;^^fe; fc) a G (0, t:) 



AE,^,^2e{n) = V5Y[{^) 



n—l 

as in (|A3|). When N is odd, though, the ground state is at the bottom of a one-paramenter, 
continuous, gapless (in the N ^ oo limit) band of eigenvalues ||l^. Within this band 

AE{a) = e{a) - e{n) a e (0, tt) (A6) 

This may be interpreted by saying that the ground state at N odd is not the "vacuum" , but it 
contains one particle (kink) which can take up a whole band of dynamical states. 

Surface energies for open antiferromagnetic chains differ, in general, when the limit iV — > cx) is 
taken for N even (e^*'+^) or N odd (e'-"'^'). They are defined by 
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Eo{N) = eoN + eS'''^^ + o{l) {N ~* oo) (A7) 



and, for the case at hand, they can be computed exactly from (A5). For the biquadratic chain, 
one finds ||,^ 

q — q 



g(s,+) ^ 1 + 2\/5 y ^ ~ 1.6550092 (A8) 

71—1 



APPENDIX B: GROUND STATE PROOF FOR PHASE < oo > 

It has to be proven that, within the boundaries of region (|2^), the ladder ground state indeed 
belongs to the sector iV^^^ = N, therefore it is just the ground state of (||) at = — 7r/2. Here 
it will be shown that it is the lowest energy state within the class which encompasses the sectors 
{iVj^^ivj"^}^^! where and {N^^^'-^i arc arbitrarily large but kept finite as ^ c». 

In other words, only one triplet fragment, Nn"^ has a diverging size in the thermodynamic limit, 
as dictated by the analysis of section 4. It will be assumed that infj>2 so(j) = so(2)- 

From (|l|) 

n-l 

AE{N) = 1) + 4^\Ni'^;l) + f^sN^'^ - 4n/.2 - {Ei^HN; 1) + f^^N) 

i=i 

When N ^ oo, setting AE = hniAr^oo AE{N) 
Since N^'^ = TV - Ar(") - EJJi ^j'^ and e^^) > e(^'+) 

n-l 

AE>J2 {4^\n^'^; 1) - eoA^j'^) - (eo + fis)N^°^ + e^''+^ - Anfi2 (Bl) 



3 = 1 



The bracketed term is so{N^^''). Actually so(j) has been defined for j > 2, whereas in (Bl) N. 



r(i) 



3 



can be 1 for some j. The definition of so{j), though, makes sense for j — 1, too. In that case 
so(l) = ~eo > so(2), so it is also true that infj>i so(j) — so(2). Then 

AE>{n^ l)so(2) + e('''+) - (cq + fi3)N^°^ - An^i2 > nso(2) - (eo + fi3)N^°'> - 4nn2 

because e'^'''+) > so{2). Next, notice that 

TV(°) > n 
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because to create a new fragment at least one singlet must be added. So, setting N'-^^ ~ n + AN'^^^ 

AE > -(eo + M3)A7V(o) _^ ^(^^(2) - eo - ^3 - 4/i2) 
which is positive due to inequalities (p4[). 
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TABLE I. Ground state energies of the biquadratic chain 



N 


e^/\n- 1) 


1 





2 


-4 


3 


-6 


4 


-9.56155281 


5 


-11.76372382 


6 


-15.14325669 


7 


-17.45402100 


8 


-20.73101075 


9 


-23.11096357 


10 


-26.32131467 


11 


-28.74967412 


12 


-31.91290040 


13 


-34.37723072 


14 


-37.50520573 


15 


-39.99742180 


16 


-43.09794702 


17 


-45.61247178 


18 


-48.69096768 


19 


-51.22377887 


20 


-54.28417520 


21 


-56.83226912 


22 


-59.87751203 


23 


-62.43858209 


24 


-65.47094086 


25 


-68.04317437 


26 


-71.06443665 


27 


-73.64638147 


28 


-76.65798212 


29 


-79.24845576 


30 


-82.25156507 


31 


-84.84959092 


32 


-87.84517669 


33 


-90.44993829 


34 


-93.43881050 


35 


-96.04961790 


36 


-99.03246169 


37 


-101.64872625 


38 


-104.62612662 


39 


-107.24734191 


40 


-110.21980250 


41 


-112.84552974 


42 


-115.81348717 


43 


-118.44334259 


44 


-121.40717895 


45 


-124.04082624 


46 


-127.00087653 


47 


-129.63801837 


48 


-132.59457884 


49 


-135.23495129 


50 


-138.18828505 


51 


-140.83165264 



TABLE II. Sequences so(n), si(n), S2(n) up to n = 51 



n 


so(n) 


si(n) 


S2(n) 


1 


n.d 


n.d. 


n.d. 


2 


1.59372686 


8 


n.d. 


3 


2.39059029 


6 


2 


4 


1.62590091 


6.37436854 


4.34232922 


5 


2.22059333 


5.88186191 


3.76372382 


6 


1.63792389 


6.05730268 


4.35744252 


7 


2.12402301 


5.81800700 


4.07241260 


8 


1.64389669 


5.92314593 


4.36550538 


9 


2.06080730 


5.77774089 


4.19041296 


10 


1.64731963 


5.84918104 


4.37049300 


11 


2.01582361 


5.74993482 


4.24989137 


12 


1.64946076 


5.80234553 


4.37387012 


13 


1.98199387 


5.72953845 


4.28469929 


14 


1.65088229 


5.77003165 


4.37630143 


15 


1.95552965 


5.71391740 


4.30709734 


16 


1.65186786 


5.74639294 


4.37813150 


17 


1.93420653 


5.70155897 


4.32249436 


18 


1.65257406 


5.72834914 


4.37955644 


19 


1.91662630 


5.69153099 


4.33360804 


20 


1.65309340 


5.71412371 


4.38069587 


21 


1.90186291 


5.68322691 


4.34193723 


22 


1.65348343 


5.70262019 


4.38162680 


23 


1.88927680 


5.67623474 


4.34836887 


24 


1.65378146 


5.69312529 


4.38240103 


25 


1.87841138 


5.67026453 


4.35345753 


26 


1.65401253 


5.68515493 


4.38305458 


27 


1.86893114 


5.66510627 


4.35756578 


28 


1.65419392 


5.67836905 


4.38361332 


29 


1.86058371 


5.66060398 


4.36093953 


30 


1.65433783 


5.67252173 


4.38409626 


31 


1.85317541 


5.65663939 


4.36375078 


32 


1.65445307 


5.66743075 


4.38451767 


33 


1.84655490 


5.65312114 


4.36612306 


34 


1.65454612 


5.66295821 


4.38488848 


35 


1.84060215 


5.64997752 


4.36814708 


36 


1.65462179 


5.65899781 


4.38521721 


37 


1.83522066 


5.64715146 


4.36989082 


38 


1.65468372 


5.65546630 


4.38551055 


39 


1.83033186 


5.64459694 


4.37140610 


40 


1.65473470 


5.65229756 


4.38577388 


41 


1.82587089 


5.64227649 


4.37273306 


42 


1.65477689 


5.64943840 


4.38601154 
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44 


1.65481197 


5.64684553 


4.38622707 


45 


1.81802811 


5.63821937 


4.37494137 


46 


1.65484125 


5.64448340 


4.38642340 


47 


1.81456284 


5.63643558 


4.37586789 


48 


1.65486580 


5.64232250 


4.38660297 


49 


1.81135678 


5.63478964 


4.37669902 


50 


1.65488645 


5.64033817 


4.38676782 


51 


1.80838229 


5.63326611 


4.37744812 



